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We discuss recent solutions for SU(2) calorons with non-trivial holonomy at higher charge, both through analytic 
means and using cooling, as well as extensive lattice studies for SU(3). 



1. Introduction 

We first briefly summarize the essential fea- 
tures of instantons at finite temperature T, so- 
called calorons A non-trivial Polyakov loop 
at spatial infinity (the so-called holonomy), which 
in some sense plays the role of a Higgs field, re- 
veals the monopole constituent nature of these 
calorons |2I3I4| . Trivial holonomy, i.e. with val- 
ues in the center of the gauge group, is typical for 
the deconfined phase. Non-trivial holonomy is 
therefore expected to play a role in the confined 
phase (r < Tc) where the trace of the Polyakov 
loop fluctuates around zero. The properties of 
instantons are therefore coupled to the order pa- 
rameter for the deconfining phase transition. 

Finite action requires the Polyakov loop at spa- 
tial infinity to be constant which for SU(7i) can 
be parametrized by 



exp(27ridiag(/ii, ^2, ■ • ■ , ^^n))g, 



(1) 



where J2m=i Z^™ ^ ^ 9 chosen to bring Voo 
to its diagonal form, with the n eigenvalues being 
ordered as < /i2 < . . . < < Mn-i-i = 1 + Mi- 
Caloron solutions are such that the total magnetic 
charge vanishes. A single caloron with topolog- 
ical charge one contains n—1 monopoles with a 
unit magnetic charge in the j-th U(l) subgroup, 
which are compensated by the n-th monopole 
of so-called type (1,1,. ..,1), having a magnetic 



charge in each of these subgroups At topo- 
logical charge Q there are \Q\n constituents, \Q\ 
monopoles of each of the n types. The monopole 
of type m has a mass Sir'^h'm/b and its core has a 
size b/ (47rfm), with i/^n = Mm-i-i— /^m- The sum rule 
J2m=i'^rn = 1 guarantees the correct action Stt^ IQj 
for calorons with topological charge Q. We can 
use the classical scale invariance to set b—1/ kT, 
the period in the imaginary time direction, to 1. 

Perturbative fluctuations give rise to an effec- 
tive potential as a function of the background 
Polyakov loop, whose minima occur where the 
latter takes values in the center of the gauge 
group 0. Recently the non-perturbative contri- 
bution due to the calorons at non-trivial holon- 
omy was calculated [7], and when added to the 
perturbative contribution, the minima at trivial 
holonomy turn unstable for decreasing tempera- 
ture, right around the expected value of Tc. This 
lends some support to monopole constituents be- 
ing the relevant degrees of freedom which drive 
the transition from a phase in which the center 
symmetry is broken at high temperatures to one 
in which the center symmetry is restored at low 
temperatures. 

The charge one SU(n) action density, S{x) = 
— ^TiF^j^{x), has a particularly simple form [S], 

TTFl{x)=dldl log [itr(A„ • • • ^i) - cos(2^i)] , 
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(2) 
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with r,„ — \x — y,n\ the distance to the m}^ con- 
stituent monopole at c„i = cosh(27rz^„ir„i), 
= sinh(27rz^„r„), r„+i ee ri and yn+i = yi. 
The chiral fermion zero mode is locahzed |8I9| 
to the m}^ constituent monopole, provided one 
uses fim < z < ^m+i for the boundary condi- 
tion f) = e-Kp{-2Triz)iiz{t,x). "Cychng" 
through the values of z gives the distinct signa- 
ture of "jumping" zero modes through which one 
can identify well-dissociated calorons, as observed 
in lattice studies for SU(2) (TUI and SU(3) TT. 

At any topological charge Q, well-separated 
constituents can be shown to act as point sources 
for the so-called far field, that is at large dis- 
tance from any of the cores, where the gauge 
field is abelian ■ When constituents of comple- 
mentary charge [n constituents of different type) 
come together, the action density no longer devi- 
ates significantly from that of a standard instan- 
ton. Its scale parameter p is related to the con- 
stituent separation d through np'^ /b — d. A typi- 
cal example for a charge one SU(2) caloron with 
far and nearby constituents is shown in Fig. ^ 
When p <C 6 no difference would be seen with 
the action density of the Harrington- Shepard so- 
lution 1 , the gauge field is nevertheless vastly 
different, as follows from the fact that within the 
confines of the peak there are n locations where 
two of the eigenvalues of the Polyakov loop coin- 
cide |13I14I15] . When, on the other hand, con- 
stituents of equal charge come together (which 
requires \Q\ > 1) an extended core structure ap- 




Figure 1. Two charge one SU(2) caloron profiles 
at t = with b = 1 and p,2 = —pi = 0.125, for 
p = 1.6 (bottom) and 0.8 (top) on equal loga- 
rithmic scales (action density cutoff at l/(2e^)). 



pears ^1. For two coinciding constituents this 
gives the typical doughnut structure also observed 
for monopoles 16 , see Fig. |21 




Figure 2. The action density (on a linear scale) of 
a typical SU(2) caloron with topological charge 2 
and p2 = —pi =0.25 for which two constituents of 
equal magnetic charge are closer than their indi- 
vidual sizes (but not exactly on top). The other 
two constituents are far away, for which this be- 
comes a static charge two monopole solution. 

2. Analytic SU(2) results at higher charge 

Also at higher topological charge the zero 
modes, denoted by ^^"(a;), where a — 1, . . . ,\Q\, 
play an important role. We can write 



n{x)^i>'i{x) ^ ~{2n)-^dlff{z,z). 



(3) 



where f°'''{z, z') is a Green's function that appears 
in the construction to be discussed below. 

2.1. The far field limit 

The trace, i.e. the sum over the zero-mode 
densities, has a remarkably simple form in the 
far field limit (denoted by ff and defined by ne- 
glecting terms that decay exponentially with the 
distance to any of the constituent cores) 



tr/*(2:, z) = 47r^V™(x), for pm < z < p,-, 



(4) 



As is implicit in the notation, Vm is static and in- 
dependent of z within its interval of definition. In 
addition Vm has to be harmonic (up to singulari- 
ties) , because the zero modes decay exponentially 
as long as z ^ pj (for any j), and therefore do 
not survive in the far field limit. For Q = 1 one 
simply has Vmix) = l/(47r|x — ym\), whereas for 
Q = 2 we found ^ 



1 



2n\x\ 



ry(if>)[ 



.(5) 
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where y{(p) ~ ( vT~-l?cos v?, 0, sincp), up to an 
arbitrary coordinate shift and rotation. Here T) 
is a scale and k a shape parameter to characterize 
arbitrary SU(2) charge 2 solutions. In this repre- 
sentation it is clear that Vm(x) is harmonic ev- 
erywhere except on a disk bounded by an ellipse 
with minor axes 2VVl^ and major axes 21?. 
Although not directly obvious, when k ^ 1 the 
support of the singularity structure is on just two 
points, separated by a distance 21?. Taking an 
arbitrary test function /(.t) one proves that |T2] 



hm / f{x)djV„,{x)<fx 



E /(0>0,^). (6) 



z=±V 



Monopoles of different charges have to adjust 
to each other to form an exact caloron solution, 
such that k and T) are in general not independent. 
So far we constructed two classes of solutions il- 
lustrated in Fig. 13 for both of which a large value 
of T) implies that k approaches 1 exponentially. 
Hence we find point-like constituents, a necessary 
requirement to describe the field configurations at 
larger distances in terms of these objects. When 
all constituents of other types are sent to infinity, 
we recover the exact multi-monopole solutions of 
a given type of magnetic charge, and our results 
therefore also provide explicit solutions for the 
monopole zero modes, which were not known be- 
fore for the multi-monopole configurations. 

y 
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Figure 3. Location of the disk singularities (shad- 
ing according to magnetic charge) for a so-called 
"rectangular" (left) and "crossed" (right) config- 
uration, as used for the Q = 2 solutions shown in 
resp. Fig. Inland Fig.^ Curves for the "crossed" 
give a one parameter family of solutions inter- 
polating between two axially symmetric solutions 
for which k = 1 independent of V (see Ref. |12j). 



Surprisingly, the charge distribution giving rise 
to the abelian field far from any of the constituent 
cores (even when extended due to overlap) can be 
calculated exactly from Vm{x)- For SU(2) we can 
parametrize the holonomy by Voo = exp(27r«a; • r) 
(^2 = —^J■l = Ml), where Tq are the usual Pauli ma- 
trices. The field asymptotically becomes abelian, 
and is necessarily proportional to w • t. Hence 
Aq{x) — 2TTiuj ■ f ~ ^iuj ■ f^{x), where the con- 
stant term shows again how Aq plays the role of a 
Higgs field. For Q = 1 and 2 we have shown that 
$(f) = 27r(Vi(f) - V2(^)). Therefore, the singu- 
larity structure in the zero-mode density agrees 
exactly with the abelian charge distribution, as 
given by df^{x). Outside the cores there is no 
source for the abelian field, implying Aq to be 
harmonic in the far field (as is also clear from the 
relation to V„i), and we find for the action density 



S^{x) = -2Tr(a,Ao(f))2 = -df^^{x). 



(7) 



Our expectation is that the result for Aq{x) will 
hold for arbitrary Q and will generalize to SU(n) 
with each Vm{x) simply associated to the U(l) 
generator with respect to which the magnetic 
charge of the monopole of type m is defined. 

2.2. The construction in brief 

A charge Q caloron with non-trivial holonomy 
is found with the Atiyah-Drinfcld-Hitchin-Manin 
formalism |17| by placing \Q\ instantons in the 
time interval [0,6], which are repeated for each 
shift of t over b, after a color rotation with Vcxi. 
The resulting infinite topological charge is dealt 
with by Fourier transformation, relating it to the 
Nahm transformation for calorons [2]- The vari- 
able z we introduced in formulating the general- 
ized boundary conditions for the chiral fermions 
is precisely the dual of t under this transforma- 
tion. Singularities are introduced at z = 
through the powers of Voo , as is seen from writing 
Voo = J2m, e^^^^^'Pm in terms of the n projectors 
P„i. The Fourier transformation also produces 
A{z) as a U(|Q|) gauge field on a circle, which 
satisfies the Nahm equation 

4-,Aj{z) + [Ao{z),Aj{z)] + ie,fe4ifc(z),i,(z)] 

= 27rjEm^(^-M™)A^, (8) 
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where 27^ClP™C6 = l2S'™-T•p^^ with Cb two-com- 
ponent spinors in the h representation of SU(n). 
Between the singularities A{z) is constant when 
Q — 1, and for (3 = 2 it can be given ^2] in terms 
of the Jacobi eUiptic functions suk, cuk and duk- 
The Green's function fx[z,z') (cmp. Eq. (Q) 
is given by g'^{z)fx{z, z')g{z'), defined through 

+ ^(^; ^)}Mz, z ) = A^H{z-z')l, (9) 

withy(z;x) = 47r2i?2(2;;f)+27r^^5(z-Ai„0S'™, 
S') = — {2TTi)^^ g{z)Aj{z)g\z) and S*™ = 
g{Hm)Smg'^ [fJ^m) ■ Notc the latter play the role of 
"impurities". The "dual" gauge function g{z) = 
exp(27ri(^o — a^ol)^), defined in terms of the dual 
holonomy exp(27ri^o) = Pexp(Jg ^o(z)d2;), can 
be used to transform Aq — 27rixol to zero, in or- 
der to simplify as much as possible the Green's 
function equation. Although fx{z,z') is periodic 
in z and z' with period 6, fx{z, z') no longer is. 

Given a solution for the Green's function, there 
are straightforward expressions for the gauge 
field mi only involving the Green's function eval- 
uated at the "impurity" locations and for the 
fermion zero modes |12| . For the zero-mode den- 
sity see Eq. As an example we give the 
Green's function at z' — z, which formally can 
be expressed as (the x dependence of Tz is im- 
plicit) fx{z,z) = -47r2 {{l~Tz)-^)^^, where the 
(1,2) component on the right-hand side of the 
first identity is with respect to the 2x2 block 
matrix structure, each of size |Q| x |Q|, and 



with T„j, Wm{z) and given by 



(■z + l 

J^2=.gT(l)Pexp / dw 



1 

V{w;x) 



(10) 



This has allowed us to find a compact formula for 
the action density jl9| . 

Tri^2^(x) = a^a^logdet {ie-^'^°{l^Tz)), (11) 

which is actually independent of z and reduces 
for Q=l to Eq. ©. 

The formal expression for Tz can be simplified 

with Tz = Wm {z)W-\ii^)T^,^ Wm {pL,n)W-\z) , 

and decomposing Tf^^ into "impurity" contribu- 
tions at z — /i„j and the "propagation" between 
them, resp. T„i and i?™ = VFm(/im+i)PF"^(^,„), 



J. 77 



'Aij' "^^^'^-[^f+iz) i-zf^(z))' 

(12) 



dzfmiz), 

•^Mm ■ ■ ■ HiTig'^{l)Hn ■ ■ ■ Tm+lHjm 

The matrices f^{z), defined for z G (/Zm, /im-i-i), 
satisfy f^{z) exp (±27r|f|(z - ^„)) /„(;«„) 
for \x\ oo. Together these form the 2\Q\ solu- 
tions of the homogeneous Green's function equa- 
tion. To resolve the full structure of the cores 
for Q = 2 we had to find the exact solutions for 
fm{z) on each interval. For this task we could 
make convenient use of an existing analytic result 
for charge 2 monopoles [21], adapting it to the 
case of calorons. Essential is that once f^{z) are 
known, everything else can be easily determined 
in terms of these, cmp. Eq. H12|) . A sample of the 
results m] can be foimd in Fig. |2 and Fig. 0] 




Figure 4. In the middle is shown the action den- 
sity in the plane of the constituents at i = for 
an SU(2) charge 2 caloron with TyVoo = in 
the "crossed" configuration of Fig. On a scale 
enhanced by a factor IOtt^ are shown the densi- 
ties for the two zero modes, using either periodic 
(left) or anti-periodic (right) boundary conditions 
in the time direction. 

3. Cooling studies for SU(2) 

Dissociated calorons in a dynamical setting 
were first observed for SU(2), using cooling to 
analyze Monte Carlo generated lattice configura- 
tions just below Tc Yet, applying the same 
methods for a symmetric box, corresponding to 
low temperatures, no dissociation was seen |15|. 
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in contrast to what was suggested by a fermionic 
zero- mode study Motivated by this we re- 

cently performed lattice cooling studies for SU(2), 
both at finite and low temperatures, to analyze 
the constituent nature of self-dual configurations 
in more detail ^22^ . At finite temperature this can 
be compared with the analytic results, but such 
results are not available for the symmetric torus. 
We used so-called e-cooling i.e. cooling with 

the improved action ( •-^ = U^{x) £ SU(n)) 



S{e) 



■ Rc Tr 1 



48 



■ Re Tr 1 - 



(13) 



The O(a^) lattice artifacts are proportional to e, 
and cause instantons to shrink for the Wilson ac- 
tion, which is recovered by putting e 1. This 
shrinking can be stabilized with e = or re- 
versed with £ < 0, which was therefore called 
over-improvement |23j. The are introduced 
to implement finite temperature (with = ^^^^ 
and = on a symmetric lattice, using an 

anisotropic coupling ^. In this way the temper- 
ature can be lowered in arbitrarily small steps, 
after each of which we perform some cooling ad- 
justing the configuration to a solution at this tem- 
perature. We called this process "adiabatic" cool- 
ing. It allows us to start from a finite temperature 
configuration with well localized constituents, in- 
vestigating how these evolve when lowering the 
temperature. 

3.1. From finite to "zero" temperature 

Constituents called instanton quarks were con- 
jectured to play a role long ago j21]. Based on 
our experience with calorons, it is expected that 
when lowering the temperature the constituents 
will grow in size. In a finite volume, as used 
on the lattice, the maximal constituent separa- 
tion is limited. We conclude that when b ^ L, 
i.e. for the symmetric box (which will be denoted 
by "zero" temperature), constituents necessarily 
overlap. In case of maximally non-trivial holon- 
omy where all types of constituents are of equal 
mass, a periodic array of constituents is actually 



related to twisted instantons of fractional topo- 
logical charge 25'T6] . As soon as the constituents 
are not equidistant, they tend to merge and can 
no longer be identified from their action density 
profiles. In Fig. [S] we present an example start- 
ing with two well-separated constituents at finite 
temperature (using ^ = 4 on a IB** lattice), low- 
ering the temperature by reducing ^ to 1. 

There is one minor complication, that on a 
torus no charge one self-dual solutions with finite 
size can exist, as can be proven with the help of 
the Nahm transformation This means that 
charge 1 instantons shrink (equivalent to the con- 
stituents getting closer) even when cooling with 
an action that has no lattice artifacts. This is why 
over-improvement was used in Fig. [S] to counter- 
act this obstruction. In Figs.Elanddthe result for 
the adiabatic cooling of a charge 2 caloron which 
is free from such obstructions is shown. For Fig. 
we start at finite temperature with the so-called 
"crossed" configuration (see Fig. 13 right). 




Figure 5. Starting from a continuum caloron so- 
lution with well-separated lumps, discretizcd on 
the anisotropic lattice and adjusted to the finite 
volume by 100 e = — 10 cooling sweeps, we re- 
duced ^ from 4 to 1, through ^ = 2i/2, 2 and \/2, 
applying between each of the 4 steps 100 e = —10 
cooling sweeps. Shown is the action density (top) 
and the the Polyakov loop in the time direction 
(bottom) along a line through the constituent lo- 
cations. The dotted, dashed and full curves are 
for C = 4, 2 and 1. 
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In Fig. the starting point at finite tempera- 
ture (top) is the "rectangular" configuration (see 
Fig. O left) . It has Tr Voo = and was taken 
from a Monte Carlo generated configuration at 
f} = 2.2, first cooling down to slightly above the 
2 instanton action with e — 1. After this, many 
thousands of e = —2 cooling sweeps were per- 
formed, followed by 500 with e = 0. We only 
show the action density in the plane through one 
of the two doughnut structures, the other dough- 
nut is separated by half a period along the axis 
of symmetry. (The periodic zero modes are local- 
ized to one doughnut, and the anti-periodic zero 
modes to the other jl2l22| '). To reach this con- 
figuration, we made use of the fact that under 
over-improved cooling constituents with opposite 
magnetic charge "repel" whereas those with the 
same magnetic charge "attract" , giving rise to 
constituents of equal magnetic charge to be ex- 
actly on top of each other. 

The adiabatic cooling is performed with e = 
cooling to prevent such "forces" due to the lat- 
tice artifacts while lowering the temperature. But 
when applying cooling with e — —10, the "fat" 
doughnut (see Fig. [7| right) ultimately will reach 
the self-dual constant curvature solution that is 
allowed ^28 29 for a symmetric box with Q = 2. 

There is another way to avoid the Q ~ 1 ob- 
struction on a torus, namely by using twisted 
boundary conditions |25| . as was used for the 
study of calorons before |13| . Under much pro- 
longed over-improved cooling with a twist in the 
0-i planes we found somewhat to our surprise that 
it was possible to push the constituents further 
apart than half a period, at the same time invert- 
ing the magnetic charge of one constituent with 
respect to the other (this is in part due to the fact 
that the Polyakov loop is anti-periodic in the i-th 
direction due to the twist, for more details see 
Ref. 122] )• This ultimately gave a perfect dough- 
nut, i.e. two magnetic monopoles of equal charge 
on top of each other. 

All these studies have shown that at low 
temperatures constituents are not localized to 
better than half the size of the volume, both for 
the action density and chiral zero modes. We 
conclude that well-localized lumps at zero tem- 
perature for low-charge self-dual backgrounds can 



only be found as instantons, even though it is 
clear that these are built from constituents of 
fractional topological charge, as seen through the 
behavior of the Polyakov loop. A good measure 
for the zero-mode localization is the inverse 
participation ratio, / = V^^p^ix), where p{x) 
is the zero-mode density '^\{x)'^ z{x) . The bigger 
/ is, the more localized is the zero mode. This can 
be contrasted with a constant zero-mode density 
for which 1 = 1. On average, at finite tempera- 
ture |TT] I is indeed considerably larger than at 
zero temperature |21) . but in the latter case val- 
ues of / as big as 20 or more are still seen to occur 
for cases where zero modes jump over distances as 
large as half the size of the vohime when cycling 




Figure 6. Starting at finite temperature with a 
charge 2 caloron in the "crossed" configuration on 
a lattice with ^ = 4 after 1000 cooling sweeps 
with £ = (left two figures), changing ^ through 
2\/2, 2, a/2 to reach 1 (right two figures), at each 
of these applying 1000 coofing sweeps with e = 0. 
The top plots show the action density integrated 
over time, the bottom plots the Polyakov loop in 
the time direction, both in the y-z plane at a; = 8 
(where all constituents lie). 




Figure 7. As in Fig. El but now starting from 
a "rectangular" configuration first made at finite 
temperature (^ = 4, left) as described in the text. 
The action densities (^ = 1, right) integrated over 
time, is shown in the plane of the doughnut. 
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through the boundary conditions. For our low- 
charge self-dual backgrounds the zero modes asso- 
ciated with separated constituents never reached 
values of / above 2. 

The possibility that the zero modes are local- 
ized to instantons (formed from closely bound 
constituents) and jump between well-separated 
instantons, rather than well-separated isolated 
constituents, was discussed in Ref. [30] . It was 
found that at finite temperature this is unlikely 
to occur, but it could not be ruled out for zero 
temperature. In a random medium of topological 
lumps the mechanism of localization of the zero 
modes could very well be similar to Anderson lo- 
calization. In such a case one perhaps should ex- 
pect a dependence on the fermion boundary con- 
ditions, even when constituents remain well hid- 
den inside instantons. In the case that instantons 
form a dense ensemble, this is similar to the state- 
ment that it is impossible to determine which set 
of constituents forms an instanton. That some 
of the zero modes, if associated to fractionally 
charged lumps, are more localized than we ob- 
served 22 can have a dynamical origin. 

3.2. Cooling histories 

That for temperatures just below Tc well disso- 
ciated calorons are seen from the action density 
profiles has now been well established ^U] • Never- 
theless, it is interesting that this can be deduced 
also from inspecting the cooling histories obtained 
with the Wilson action. There are three processes 
under which the action is lowered during cooling. 




The first and most important one is the removal 
of short distance fluctuations. Next is the annihi- 
lation of lumps with opposite topological charge, 
and last the falling of instantons through the lat- 
tice due to its shrinking. Annihilation can, how- 
ever, also occur between well-localized lumps of 
opposite fractional topological charge. 

To detect these events one defines a plateau 
in the cooling history as the point of inflection 
for the action as a function of the cooling sweeps 
(i.e. when the decrease per step becomes mini- 
mal) jlOj . In the confined phase, where on average 
the charge fraction is one half, the change in ac- 
tion would typically be around one unit, and the 
topological charge would remain unchanged (an 
example is indicated by "s" in Fig. ^ . This can 
thus be easily distinguished from the annihilation 
of an instanton and anti-instanton, for which the 
action changes by two units. On the other hand, 




Figure 8. Sample of 8 histories on a 16"^ x 4 lattice 
at /3 = 2.2 (T « O.STc). Pluses give the Wilson 
action and crosses the absolute value of the (order 
improved clover averaged) topological charge. 



Figure 9. Scatterplots of |A(3| (horizontally) ver- 
sus AS* (vertically) for a G'^ x 4 lattice at (i) 
/3 = 2.4 (T w 1.2rc), (ii) /3 = 2.2 (T w CST^), 
and (iii) a 16^^ lattice bX 13 ^ 2.3 (T « 0.25rc). 

when an instanton falls through the lattice (i.e. 
two constituents with the same sign for their frac- 
tional topological charge, but of opposite mag- 
netic and electric charge, come together and form 
a localized instanton, which then shrinks under 
cooling) both the topological charge and the ac- 
tion changes by one unit. 

We collected the information concerning cool- 
ing histories in scatterplots, each based on 50 
Monte Carlo generated configuration, see Fig. |5| 
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Only for T just below the constituents are suf- 
ficiently well localized to give rise to clear con- 
stituent annihilations as seen from the clustering 
of (AS, lAQI) around (1,0), measured from one 
plateau to the next. 

The presence of constituents of opposite topo- 
logical charge also gives rise to the interesting fact 
that the number of near zero modes (in a smooth, 
but non-selfdual background) can depend on the 
fermion boundary conditions as shown in Fig. 1101 
The anti-periodic zero modes localize to con- 
stituents for which the Polyakov loop is — 1, 
whereas for periodic zero modes P — The 
underlying configuration of Fig. IIUI has 3 anti- 
selfdual lumps of which one has P = -1-1 and one 
self-dual lump with P = — 1 . Thus there are three 
anti-periodic near zero modes, two of which dis- 
appear as soon as z ^ \. 



cycling be zero mode + 
periodic be non-zero modes x 
cycling be near-zero modes Q 
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Figure 10. Zero-modes for an anti-caloron on a 
\& X 4 lattice, with in addition two constituents 
of opposite fractional topological charge. Crosses 
give the low-lying eigenvalues A of the clover- 
improved Wilson-Dirac operator with periodic 
boundary conditions in time. The curves trace 
two near zero modes from their value with anti- 
periodic (left) boundary conditions to the peri- 
odic (right) case, in 5 equal steps in z (squares). 

4. Marginally stable configurations 

The last cooling run shown in Fig. |H1 ends in 
a configuration with (5 = and S" w 1/4 units, 
which is stable under further cooling. These so- 
called Dirac sheets had been shown to typically 



arise after constituent annihilation, leaving a con- 
stant magnetic field |81j . It is perhaps surprising 
that such constant magnetic fields can be sta- 
ble, but this is a consequence of the finite vol- 
ume in which a non-trivial holonomy can affect 
the fluctuation spectrum. Since the action of the 
constant magnetic field does not depend on the 
holonomy, which also holds on the lattice, this 
was called marginal stability \62\ and explains the 
stability under further cooling, once it is reached. 

The finite temperature setting is essential, as 
otherwise one would need to introduce twisted 
boundary conditions to construct marginally sta- 
ble configurations [32]. For SU{2) the constant 
magnetic field solutions are in fact abelian, with 
= T^iT^n^^KL^L^) in a suitable gauge. Here 
rifj^i, is antisymmetric, taking even integer values 
(in the absence of twisted boundary conditions). 
Adding a constant to the abelian vector potential 
can be partly, but not completely, absorbed by 
a translation. The part that cannot be absorbed 
will affect the fluctuation spectrum. With Li =L, 
noi, — and rrii = \£ijknjk, the eigenvalues for the 
quadratic fluctuations can be shown to be ;29»33] 
\±=2Ti{2n+l±2)\rh\/ L\{2'KepMh-Cf/ {\m\Lf+ 
{2'Kq+CaY /b'^ and Aq = {2iTk^/L^Y, with multi- 
plicities resp. 2|e| and 2, where e is the greatest 
common divisor of the m;. All quantum numbers 
(n,p, g, k^) are integer (but n > 0). 

At finite temperature the spatial holonomies 
are typically trivial and the relevant Polyakov 
loop to consider for the dependence on the spec- 
trum is Po = iTrexp(iCoT3/2) = cos(Co/2) or 
rather Pq. This is because the fluctuations in- 
volve adjoint fields and the spectrum is periodic 
under a shift of Cq over 27r, whereas Pq is anti- 
periodic. To find the lowest eigenvalue we put 
n = q=pi~C'i~0 (we may restrict |Co| < tt) such 
that A_ — — 27r|?7i|/L^ -I- Cq/IP' is negative unless 
the Polyakov loop Pq is sufficiently non-trivial. 
All eigenvalues are therefore positive when 

H<^ and |Po|<cos(^). (14) 

These conditions can never be satisfied when 
b = L, nor in the deconfined phase where the 
Polyakov loop is too close to being trivial. We 
found in the cooling histories on 16"^ x 4 lat- 
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tices stable configurations witli one, two and tliree 
quarters of the instanton action (with one, two or 
three of the components of m being 2). For the 
minimal case, \rh\ — 2, the correlation between 
stability and the value of the Polyakov loop was 
studied in much detail , showing perfect agree- 
ment with the analytic prediction in Eq. (|14|l . 

5. Lattice calorons for SU(3) 

Caloron solutions with non-trivial holonomy 
for SU(3) are qualitatively different from those 
with trivial holonomy and from embedded SU(2) 
solutions, which appear only as limiting cases. 
We will illustrate this with generic SU(3) solu- 
tions obtained by cooling on a finite lattice. 

First we consider a Q = 2 solution on a 20'^ x 4 
lattice, see Fig. ^2 Eigenvalues of the Polyakov 
loop all agree very well in the region where the 
action density is lowest (10 % of the volume), and 
are used to define the holonomy Voo in the case 
of a finite lattice. For calorons with one unit of 
topological charge it can be shown that at the 
positions where two eigenvalues of the Polyakov 
loop coincide, these eigenvalues are given resp. by 

|g-iri;i3 g-iri/^a g2iri/23 1 |g27ri;ii ^-Trifii g-7ripi| 

and {_e-'^V2^g27riM2^ _g-"P2} xhis is tested on 
the right in Fig. ^2 The typical range for P (one- 
third the trace of the Polyakov loop) is shown for 




Figure 11. Shown on the left by colored spheres 
are the six constituent monopoles for a Q = 2 
caloron. The behavior of the Polyakov loop is 
shown on the right, with the filled circle corre- 
sponding to the holonomy and the three open 
circles to the predicted values of P at the con- 
stituent locations. The actual values of P for 
points in their neighborhood are represented by 
the colored corners. 



two hypothetical Q = 1 objects. It is spanned 
by the colored corners which correspond to the 
actual values inside the cores of the constituent 
monopoles. The colors distinguish the three dif- 
ferent types of monopoles. For clarity we do not 
show the full scatterplot for P. The boundary 
for the range of P values in the complex plane 
is indicated by the curved triangle. The close- 
ness of two of the three eigenvalues is represented 
by the distance to this boundary. The six iso- 
surfaces of that distance function are pictured in 
corresponding colors on the left in Fig. ^2 this is 
also where the action or topological charge den- 
sity (not shown) is maximal. 

The two zero modes for the charge 2 configu- 
ration of Fig. ^] are shown on the left of Fig. 1121 
as isosurfaces of its zero-mode density. Using col- 
ors we combine in one figure three values of z 
that specify the fermion boundary condition in 
the time direction. These z values correspond to 
the three maxima of the inverse participation ra- 
tio 51] / (shown on the right), for which the zero 
modes are maximally localized. The monopole 
locations defined through the coinciding eigen- 
values of the Polyakov loop agree (within one 
lattice spacing) with the maxima of the zero- 
mode density. With help of the z-dependent zero 
modes one can focus on one of the three types of 
monopoles alone, which brings out more clearly 
the constituent nature of calorons. 

Compared with the analytic caloron solutions 
with arbitrary holonomy for 5^ x IR"^, lattice solu- 
tions have imperfections, some specific for SU(3). 
These are caused by the finite volume with peri- 




Figure 12. The positions (red, blue, green) of the 
two zero modes (left) for the charge 2 configu- 
ration shown in Fig. 1111 with the z values corre- 
sponding to the maxima of / (shown on the right). 
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odic boundary conditions (leading to the charge 
one obstruction discussed in Sect. ^ and its fi- 
nite aspect ratio b/L, as well as by lattice ar- 
tifacts. Hence, several compromises are neces- 
sary to describe the moduli space of caloron solu- 
tions using cooling on the lattice. In this spirit we 
have performed a systematic investigation j34j of 
the structure of SU(3) calorons with \Q\ = 1,2,3 
and 4 on x Nt lattices with Ns = 12 and 20 
{Nt > 4). For Ns = 12 the results are based 
on 0(10000) cooled configurations. The holon- 
omy is non-trivial when applying cooling to con- 
figurations generated in the confined phase us- 
ing f3 = 5.6 (corresponding to T w 0.95Tc for 
Nt = 4). Similar to what was seen for SU(2) [rUj . 
the distribution of the Polyakov loop broadens 
with cooling, and cuts in Poo = Tr Voo/i are some- 
times helpful to highlight the constituent nature. 

For an automated scanning of the solutions we 
introduced two observables, the non-staticity TU' 
6t = l'^(^) —s{x + i) \ /S and the violation 

of (anti)self-duality (5_F = |s(a;) — 87r^|9(a:)| | /S. 
The non-staticity is used to determine when sep- 
arated lumps are seen within a solution. Points 
of bifurcation occur when 51 — 0.27 and 5** = 
0.06 for which respectively two and three peaks 
become visible in the topological charge den- 
sity. These values were estimated on the basis 
of the exact charge one caloron solutions with 
maximally non-trivial holonomy and suitably ar- 
ranged constituents [SI]. The C'(a'*)-improved 
field strength "HF was used to construct the ac- 
tion and topological charge densities, to minimize 
the effect of lattice artifacts. For the fermions 
with adjustable boundary conditions we used the 
clover-improved Wilson-Dirac operator. 

Like for SU(2), instantons and calorons shrink 
under cooling with the Wilson action. In addi- 
tion, charge one caloron solutions can only be ob- 
tained at the expense of violation of (anti)sclf- 
duality. For example, cooling on a 16^ x 4 lattice 
with a stopping criterion based on a minimal vi- 
olation of the lattice equations of motion JHIj or 
a minimal value oi 5p^ one is unable to produce 
dissociated |Q| = 1 calorons, although this can be 
corrected for using over-improved cooling, as we 
have seen for SU(2) in Sect. |3| However, for our 
high-statistics study we have restricted ourselves 



to Wilson-type cooling. Instead we modified the 
condition under which cooling is stopped and the 
approximate solution is stored. 

The extended stopping criterion is based on 
either the non-staticity 5t or the violation of 
(anti)self-duality Sp to pass through a minimum. 
The distribution of (almost) classical solutions on 
a 12^ X 4 lattice with respect to 5t is compared 
for both cases in the upper row of Fig. E| The 
vertical line, drawn to guide the eye, represents 
51 where the bifurcation of 1 into 2 peaks occurs. 
With the modified stopping criterion, which al- 
lows for some violation of (anti)self-duality, many 
perfectly static calorons are found which are 
clearly split into constituents. These are marked 
by "(a)" on the top right plot in Fig. to be 
contrasted with the behavior in the top left plot, 
and are mainly associated to \Q\ = 1 (with some 
IQI = 2) configurations. Yet all satisfy Sp < 0.2, 
and are sufficiently close to a classical solution. In 
the lower row of Fig.^jthe distribution over 5t is 
shown for approximate solutions found on lattices 
12^^ X 6 and 12'* (the symmetric torus), again using 




non-staticity non-staticity 8^ 



Figure 13. Histograms of non-staticity St for 
configurations obtained when the violation of 
(anti)self-duality 5f minimal (upper left) and 
when using the extended stopping criterion de- 
scribed in the text (upper right), both on a 12'^x4 
lattice. The bottom row shows the results for 
12^^ X 6 (left) and 12^* (right) lattices, also using 
the extended stopping criterion. 
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the extended stopping criterion. This confirms for 
SU(3) the tendency towards non-dissociated con- 
stituents also found for SU(2) iSj where it is now 
well understood [221 j as was discussed in Sect. 13 
As emphasized before, constituent monopoles 
can be detected through the behavior of the 
Polyakov loop, even when they are close to- 
gether. This is exemplified by Fig. 1141 where the 
non-staticity is plotted versus the average inter- 
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Figure 14. Correlation between the non-staticity 
St and the average inter-monopole distance for 
charge one configurations on 12'^ x 4. 

monopole distance for charge one calorons on 
12'^ X 4 lattices. The two clusters marked "(a)" 
and "(b)" agree with the configurations associ- 
ated to the two peaks marked in Fig. ^| for the 
top right plot, based on the extended stopping 
criterion. 

It is interesting to see how Sp is actually dis- 
tributed. This is shown in Fig. El (where a cut in 




0.05 0.1 0.15 

violation of self duality 5p 



Figure 15. Histograms for the violation of (anti) 
self-duality (5i? , separately for IQj = 1 and |Q| > 1 
on a lattice of size 12'^ x 4. 



Sp at 0.2 was applied) contrasting |(5| 1 with 
IQI > 1. Absence of configurations with charge 
one for Sp < 0.016 illustrates the non-existence of 
self-dual fields with \Q\ = 1 on the torus. 

Regardless of the number of lumps, our method 
of locating and counting monopoles described 
above is reliable provided the configurations have 
non-trivial holonomy. The number of monopoles 
is distributed very narrowly around 3\Q\. This 
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Figure 16. Average number of monopoles and 
of maxima of the topological density and their 
standard deviations as function of \Q\ for the 
12^ X 4 lattice, without cuts (left), with a cut 
|.Poo| < 0.2 (center) and with an additional cut 
St < 0.2 (right). 

is shown in Fig. 1161 where the average number 
of monopoles and peaks in the topological charge 
density are represented as a function of \Q\, with- 
out cuts and with cuts in |Poo| and St- These cuts 
effect the number of peaks but not the number 
of monopoles, which however does become more 
narrowly distributed. 

We have thus seen that the approximate SU(3) 
caloron solutions produced by cooling, even on 
relatively small lattices, resemble the analytically 
known caloron solutions. The way they manifest 
themselves depends on the aspect ratio of the lat- 
tice. 

6. Conclusions 

We have seen that instantons at finite tem- 
perature are composed of constituent monopoles. 
More details for SU(3) will appear soon (23 ■ The 
hope is of course to be able to develop a calcula- 
tional scheme to address questions like monopole 
condensation in the popular scenario of the dual 
superconductor From the analytic side this 
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is complicated by the fact that at low temper- 
ature the coupling constant is large, and instan- 
tons form a dense ensemble. More detailed lattice 
studies would be welcome. A dense ensemble of 
instantons implies that the monopoles are dense 
as well. The confining electric phase could per- 
haps be characterized by a dual deconfining mag- 
netic phase, where the dual deconfinement is due 
to the large monopole density, similar in spirit to 
high density induced quark deconfinement. 
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